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Abstract 

A method is developed to predict the vibration 
response of a composite panel and the resulting far- 
held acoustic radiation due to acoustic excitation. The 
acoustic excitation is assumed to consist of obliquely 
incident plane waves. The panel is modeled by a finite 
element analysis and the radiated held is predicted us- 
ing Rayleigh’s integral. The approach can easily in- 
clude other effects such as shape memory alloy (SMA) 
fiber reinforcement, large deflection thermal postbuck- 
ling, and non-symmetric SMA distribution or lamina- 
tion. Transmission loss predictions for the case of an 
aluminum panel excited by a harmonic acoustic pres- 
sure are shown to compare very well with a classical 
analysis. Results for a composite panel with and with- 
out shape memory alloy reinforcement are also pre- 
sented. The preliminary results demonstrate that the 
transmission loss can be significantly increased with 
shape memory alloy reinforcement. The mechanisms 
for further transmission loss improvement are identi- 
fied and discussed. 
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Introduction 


Interior or cabin noise is one important issue in 
the development and design of advanced subsonic and 
supersonic aircraft. Conventional aircraft typically 
employ passive treatments, such as constrained layer 
damping and acoustic absorption materials, to reduce 
the structural response and resulting acoustic levels 
in the aircraft interior. These techniques require sig- 
nificant addition of mass and only attenuate relatively 
high frequency noise transmitted through the fuselage. 
Improved reduction in the transmitted noise at the 
lower frequencies is required for advanced subsonic air- 
craft. Reductions will also be imperative, along with 
the need for improved fatigue resistance without sig- 
nificant mass addition, for proposed supersonic trans- 
ports, where the excitation levels will be high due to 
engine noise and turbulent boundary layer fluctuating 
pressures. 

Although structural acoustic coupling is in gen- 
eral very important in the study of noise transmission 
through an aircraft fuselage, analysis of noise trans- 
mission through a panel supported in an infinite rigid 
baffle (separating two semi-infinite acoustic domains) 
can be useful for studies involving active /adaptive ma- 
terials, complex loading, etc. Transmission of noise 
through conventional baffled panels has been studied 
extensively 1 . More recent work has been aimed at de- 
veloping active and/or adaptive methods of controlling 
the structural acoustic response of panels to reduce 
the transmitted noise 2 . These methods have focused 
on using piezoelectric materials, and to a lesser extent 
shape memory alloys, to control the structural acoustic 
response of panels subject to harmonic acoustic exci- 
tation. 

Turner et. al. 3 developed a finitf element formu- 
lation to study the response of shape memory alloy 
hybrid composite panels due to thermal and acoustic 
loads. A finite element approach was used to allow for 
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the SMA material nonlinearity, thermal expansion ef- 
fects, geometric nonlinearities due to thermal buckling, 
and complicated geometries and boundary conditions. 

The purpose of the present work is to develop 
an approach for analyzing the structural acoustic re- 
sponse of panels due to obliquely-incident plane-wave 
acoustic excitation of general temporal nature, with 
the intent for subsequent studies involving SMA hy- 
brid composite panels. The method is validated 
through comparison with a classical approach for har- 
monic excitation 1 . Results are also presented for com- 
posite panels with and without shape memory alloy 
reinforcement . 


Shape Memory Alloy Hybrid Composites 


Shape memory alloys exhibit a characteristic 
phase transformation from martensite to austenite, 
initiating at the austenite start temperature T s and 
asymptotically ending at the austenite finish tempera- 
ture Tf . A shape memory alloy in the low temperature 
martensitic condition ( T < T s ), when plastically de- 
formed and the external stresses removed, will regain 
its original (memory) shape when heated. For exam- 
ple, strains of typically six to eight percent can be com- 
pletely recovered by heating the nickel-titanium alloys 
(Nitinol) above the austenite finish temperature. The 
transformation temperatures can be altered by chang- 
ing the composition of the alloy. In addition, when 
Nitinol is heated, the Young’s modulus increases three 
to four times and the yielding strength also increases 
approximately ten times 4 1 . 

The materials, referred to in this study as SMA 
fiber-reinforced hybrid composites, are conventional 
advanced composite materials (such as graphite-epoxy) 
that contain embedded SMA fibers having the same di- 
rection as the graphite fibers. The memory effect of 
the SMA fibers may be put to use by applying an ini- 
tial elongation to them during fabrication of the lam- 
inate. Thus, once the laminate is heated above T s , 
the SMA fibers will try to recover the initial strain 
and induce tensile inplane forces in the laminate. The 
overall stiffness of the hybrid composite panel will bi, 
increased due to: 1) the increase of the Young’s mod- 
ulus of the SMA fibers by a factor of three or four, 
and 2) the internal tensile inplane forces induced in 
the panel from the recovery of initial strains of the 
SMA fibers. Therefore, the root.-mean-square (RMS) 
maximum deflection due to acoustic pressure loading, 
will be reduced at elevated temperatures as compared 
to composite panels without embedded SMA fibers. 

Shape memory alloys have been applied as ac- 
tuators for active control of buckling of beams 8 and 


shape control of beams 9 . It is also being studied for 
use in active vibration control of beams 10,1 1 and large 
space structures 12 . Active vibration control of flexible 
linkage mechanisms using SMA fiber-reinforced com- 
posites has been investigated by Venkatesh et al. 13 . 
Acoustic transmission and radiation control through 
the use of SMA fibers in a hybrid composite was pre- 
sented by Liang et al. 14 and Anders et al. 15 . 


Finite Element Formulation 

The finite element equations governing the re- 
sponse of a composite panel, with shape memory alloy 
wire reinforcement and subject to thermal and acous- 
tic loads, have the form 3 
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[M}{w] + ([A'] + [K*] ~ [Knat] +\[N1] 

+ ^ 2 ]) {it} = {P(t) } + { p AT } - {p; } ( 2 ) 


where [M\ and [A"] are the system mass and stiff- 
ness matrices; [A"*] and [A'jvat] are the geometric 
stiffness matrices due to the SMA recovery stress a* 
(or {A*}) and thermal inplane force: vector {Aat}, 
respectively; [Al] and [A2] are the first and second 
order nonlinear stiffness matrices which depend lin- 
early and quadratically upon displacement {IF}, re- 
spectively; {.P(f)} is the acoustic excitation load vec- 
tor, {P^t} is the thermal load vector, and { P '*} is 
the SMA recovery force vector. The subscripts b 
and m denote bending and membrane components, 
respectively; and the subscripts B, Nm and NB 
indicate that the corresponding stiffness matrix is 
due to the laminate stiffness [5], membrane forces 
{A m }(= [A]{e,°„}) and {A s }(= [5 ]{k}), respectively. 
Detailed derivations of the equations of motion and 
expressions for the element matrices and load vectors 
are referred to Locke 11 ’ 1 '. 

Consider a rectangular composite panel supported 
in an infinite rigid baffle, subjected to an incident 
acoustic pressure on one side and allowed to radiate to 
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a semi-infinite acoustic field on the other. The acoustic 
excitation load vector in equation 2 has the form 



with 

{ p b(t)} e = J p(£,,ihi){F w (£,,il)} dA e (4) 

A e 

and 

= {0} (5) 

where the F w are the out-of-plane interpolation func- 
tions, (£,r)) are the coordinates in the plane of the 
plate, and the subscript e indicates an element quan- 
tity. The element contributions from equation 4 are 
assembled to form {-Pj(f)} for the system. In general, 
the pressure loading on the panel consists of the in- 
cident. pressure, the reflected pressure, and the trans- 
mitted pressure: 

=Pi(t,V,t) + Pr(t,V,t) -Pt(Cv^) ( 6 ) 

However, p(£,r),t) can be written in terms of the 
blocked pressure and the radiated pressure, as sug- 
gested by Roussos 1 . The blocked pressure is that pres- 
sure on the incident side when the panel is consid- 
ered rigid and the radiated pressure is that due to 
panel motion. The acoustic radiation problem can 
be solved separately from the forced response of the 
panel by assuming that the radiated pressure is neg- 
ligible compared to the blocked pressure. Thus, the 
load vector {-Pj(f)} is due to the blocked pressure 
{pitvd) = Pbi^V’t) = 2 pdCV’t)) and the trans- 
mitted pressure (pressure radiated to the free field) 
can be determined from the resulting panel motion, 
without the need for a fully coupled analysis. 

Three different types of analyses are required to 
determine the response of SMA hybrid composite pan- 
els to an acoustic pressure at elevated temperatures. 
They are: (1) thermal buckling analysis, (2) thermal 
post.buckling analysis, and (3) vibration analysis of flat 
or thermally buckled panels. The governing equations 
for the linear dynamic response of such a panel, as- 
suming that the stable static configuration has been 
determined, are given by 3 

[M]{w}+ [K]{W} = {P(t)} (7) 

where [A"] is the stiffness matrix for the flat, or buckled 
panel. The natural frequencies u> r and mode shapes 
{</>,,} of vibration are obtained from the eigenvalue 
problem 

uj;[M]{<Pr}= ( 8 ) 


A set. of uncoupled modal equations with reduced de- 
grees of freedom can thus be obtained from equation 
7 as 

\im + r2(r^rj{?} + = {/wi (9) 

where r = 1 , 2 ,..., m is the number of modes retained, 
the delimiters |"J indicate a. diagonal matrix, and |"7J 
is the identity matrix, indicating a. mass normalized 
modal matrix. In Equations 9, a. modal damping term 
has been included where is the modal damping ratio. 
The physical degrees of freedom can be recovered from 
the truncated modal transformation 

m 

{w} = ^2 {Ma = MM m < N i (io) 

r — 1 

where N\ is the total number of degrees of freedom. 

The modal response, solution to the modal equa- 
tions of motion (equations 9), can be determined by 
time domain analysis or frequency domain analysis. 
In either case, a. continuous or discrete approach can 
be employed. In this study, analysis is done in the fre- 
quency domain through the use of a. discrete technique's. 
This allows complicated frequency distributions to be 
handled in a. straightforward manner. 

Discrete Spectral Analysis 

A discrete spectral analysis approach was chosen 
since continuous approaches are limited to rather sim- 
ple loading cases, which are not. very realistic for pra.o 
t.ica.l applications and not. easily generated in an ex- 
perimental situation. In practice, there is always a. 
finite amount, of data. Thus, one would employ a. Fi- 
nite Fourier Transform to a. sampled time history as 
seen through a. viewing window. Furthermore, for dis- 
cretely sampled time data, a. Discrete Fourier Trans- 
form (DFT) can be used to estimate the Finite Fourier 
Transforms, etc. A particularly efficient, means of per- 
forming a. DFT is known as the Fast. Fourier Transform 
(FFT). Of greater importance is that. most, experimen- 
tal data, is processed via. the FFT. Thus, a. means of 
using experimentally measured data, has obvious ben- 
eficial implications. 

Panel Response 

A discrete spectral analysis can be employed by 
assuming that, the modal force vector, and thus the 
modal response, can be written as the sum of harmonic 
functions: 

JV 2 -1 

{/(*»= y. !/ 1' ■ ■ (id 

n = 0 
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iV 2 — 1 

{«(<)} = e ^y wnt ( is ) 

n=0 

where # 2 is the number of discrete frequencies in the 
transformation and {/„} and {</„} are the complex 
spectrum amplitudes vectors, one for each frequency 18 . 
Substitution into the modal equations of motion, equa- 
tion 9, gives 


where the \H(ui n )\ is a diagonal matrix of frequency 
response functions given by 


Hr [■ J ~n ) 


1 

W-? - ojI + i2( r u> r u> n 


(19) 


Note that from equation 12, the displacement, velocity, 
and acceleration time responses can be computed by 
an inverse FFT, assuming that the phase information 
has been preserved, from 


n 2 -i 

E ( — m T^Cr^rJ 

n = 0 

N 2 - 1 

+ K 2 J){9n}e“" t = E !/ I' ■ ' 

n=0 

. (13) 

Therefore, for any one of the discrete frequencies aj n 
the relation becomes 


n 2 -i 

{m}= E ^y Wnt 

n=0 

N 2 - 1 

{«(*)} = E ( 20 ) 

n=0 

iV 2 -l 

{g(f)} = E 

n=0 


(|~ — 1 w nJ + r*2c,Vw r w n J + \uj, t_|){g„} — {/„} (14) 

The modal forcing vector is related to i lie acoustic 
excitation load vector from the system diagonalization, 
equation 9: 

{/(*)} =MV(m or {/„} = M T {P„} (15) 

Furthermore, the nonzero elements of that vector re- 
sult from {Pb(t )}• The incident pressure is assumed to 
consist of obliquely incident plane waves, so the inci- 
dent. pressure spectral amplitude vector for an element 
at a given frequency can be derived from equation 4 
and becomes 


{P n } e = j ■2 Pn {F w (^v)}e- i ' w ’' iA ^ 1 dA e 

l e 

= { r } e Pn 


(16) 


where {T} e contains the integrated element spatial 
description of the incident acoustic pressure, p n is 
a spectral amplitude of the acoustic pressure, k = 
sin#,; cos <f>i, and A = sin 0,; sin (/>,; , see figure 1. There- 
fore, assembly of the element contributions from equa- 
tion 16 results in 


{fn} = M T {r} P „ (17) 

The relation in equation 14 is simply a set of uncou- 
pled algebraic equations for the modal displacement 
spectral amplitudes. Thus, the modal displacement 
spectral amplitude vector becomes 


If the excitation is deterministic, the panel dynamic 
response solution is complete. However, for stochastic 
excitation, a statistical approach must be used. Ad- 
ditionally, panel response cross-correlation quantities 
are of interest for the acoustic radiation discussion to 
follow in the next section. Therefore, the necessary 
relations for stochastic panel response are developed 
subsequently. 

The modal displacement response correlation ma- 
trix can be formed from equation 12 as 

[RqqtflHl)] = £[{?(fl)}{?(f2)} T ] 

JV 2 -1W 2 -1 

m = 0 n=0 

( 21 ) 

Substituting from equations 16 and 18 gives 

N 2 -1N 2 -1 

[R qq (tlj2)} = E E rtf(-"m)JMV}£[PmP»] 

m = 0 n=0 

{r 

( 22 ) 

Note that if the incident acoustic pressure is deter- 
ministic, the ensemble average indicated by the ex- 
pected value in equation 22 is not necessary. If the 
incident pressure is stochastic, the necessary ensem- 
ble averages can be performed ahead of time. For a 
stationary process, the expected value takes a special 
form 19 

R[PmPn\ — E[PmPn]$mn (23) 

Then the stationary modal displacement correlation 
matrix results 


{q n } = \H(ui n )\{f n } 


( 18 ) 


[R qq ( T )] 


n 2 - i 


n = 0 


( 24 ) 


4 

American Institute of Aeronautics and Astronautics 



The stationary modal velocity and acceleration corre- 
lation matrices follow 

JV 2 -i 

[f?,,(r)]= £ 

n=0 (25) 

JV 2 - 1 ' ' 

[ff^(r)] = £ ^{ti}{ q n} T V^ T 

n=0 

The correlation matrices with a time lag of zero 
(r = 0, so that the diagonal elements are mean square 
values) can be determined from the modal response 
power spectrum. An alternative method, which is an- 
alytically consistent (so as to avoid normalization fac- 
tors) involves performing an inverse FFT based upon 
the above equations. Then the zero time lag corre- 
lation matrices correspond to the first, element of the 
inverse transforms. Finally, the stationary response 
correlation matrices of the physical degrees of freedom 
result from the modal expansion, equation 10. 

Acoustic Radiation 


Pt(x. y. z.t) = ^ p n (x, y,z)e' w ’' t (30) 

n=0 

As in i lie;- panel response section, if the excitation is 
deterministic, an expression for the radiated pressure 
spectral amplitude could be developed and used in 
an inverse FFT to determine the radiated pressure 
time response. However, for stochastic excitation, a 
statistical approach must be used. 

Substitution of equation 29 into equation 28 re- 
sults in the following expression for the transmitted 
pressure 

Pt(x,y,z,t) = 

k=1 l 

n 2 - i 

E dA k 

n = 0 

(31) 

The transmitted acoustic pressure cross-correlation 
function follows as 


The transmitted acoustic pressure radiated by 
the panel to a point (x.y. z) in the acoustic free-held 
can be related to i lie dynamic response of the panel 
through the use of Rayleigh’s integral in the form 20 

p t (x,y,z,t) = J ^w^.y.t- 1 -^ dA (26) 

where p is the fluid density, ib is the out-of-plane ac- 
celeration response of the panel, ^ and y are spatial 
coordinates in the plane of the panel, r is the distance 
from (^, y) on the panel to (x, y. z) in the acoustic held, 
c is the acoustic wave speed, and t — is the retarded 
time between the two points. Since the physical de- 
grees of freedom are known on an elemental basis 

w le) (£,TiJ)K{F w (Z,Ti)}? e) {W(t)} le) (27) 

the area integral in the expression for the transmitted 
pressure, equation 26, is taken to be the sum of element 
area, integrals: 

p t (x,y,z,t)= -^E / 

*='L' (28) 

M' 

As before, a. discrete spectral analysis can be employed 
t.o write 

JV 2 -i 

- ;)} = E (29) 

n=0 


nel nel 


Rpp{xi,ti, X2A2) — ( 27 ]-) J J 


k = 1 1=1 
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{F w {S.^l)} T k E E ^ 


(32) 


m = 0 n=0 


e- i ^-^) e ’M^- L t) {Fw(C ^ v ' )}i dAk dA ' 


where ?q is the distance from (£,y) on the panel to 
T 1 = (a? 1 , t/i, zi) in the acoustic held and ro is 1 lit;' 
distance from (£',?/) on the panel to To = (x 2 .y 2 .~ 2 ) 
in the acoustic held. Recall that the expected value 
of the spectral amplitudes takes a. simplified form for 
stationary response: 


E[{WZ } k {W„ }f] = E[{w: n } k {W„ (33) 


Then the stationary cross-correlation function results 
Rpp(xi.ti. X2A2) — * R pp (xi.X2.t). which can be fur- 
ther simplified for zero time lag (r = 0, not. to be 
confused with zero retarded time lag). Thus, the 
mean square transmitted acoustic pressure at a. point. 
(x.y.z) becomes 


nel nel 




N- 2—1 


k = 1 l— 1 


r\ r 


1 (€, 


A k A' 


E ^tE[{W*} k {W n }f]e i ^ : ^ ±) 
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( 34 ) 
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where jq and are distances from the same point 
(x,y,z) in the acoustic field to different points on the 
panel, and (?q — r[)/c is the retarded time difference. 
The element area integrals can be conveniently evalu- 
ated via Gauss-Legendre numerical integration, which 
precludes the necessity for far-held approximations. 
However, in the discussions to follow, it will be useful 
to introduce acoustic far-held assumptions for acoustic 
intensity Calculations. These assumptions entail mak- 
ing the approximations 1 


1^1 1 _ 1 

R r[ R 


(35) 


and 




. ({'-{)K + | v'-V )A 


(36) 


where Iv = sin 9 cos <ft and A = sin 8 sin <f>, see hgure 2. 
Thus, the mean square transmitted acoustic pressure 
can be written for the far-held as 


nel 


Pt(R,o,<t>) = 


‘2irR 


EE 

k= 1 1=1 




A k A' 


N 2 - 1 


E “^[{wZhiWnify' 

n = 0 


(i'-i)K+(i'-qA 


{^(£W)}, dA k d.A\ 

(37) 

Note that the interpolation functions F w are usually 
known in terms of element local coordinates. There- 
fore, the complex exponential function in the integrand 
of equation 37 must be separated into two components: 
one which is constant for each element index pair (k, /); 
the other must be included in the integration over 
the element areas in local coordinates. The area in- 
tegrations can be performed numerically, via Gauss- 
Legendre numerical integration, but significant reduc- 
tions in computational time can be realized by per- 
forming the integrations in closed form. Other acous- 
tic quantities, such as the acoustic intensity and total 
radiated power, can be related to this mean square 
pressure in equation 37. 


The acoustic intensity at a point is equal to the 
time average of the product of the perturbed pressure 
and velocity. For a pressure and velocity that are ran- 
dom processes, this temporal average must be taken 
in the limit as the period T — *■ oo: 


I{x,y,z) = (p{x,y,z,t)u{x,y,z,t)) 


1 

= Jim l p{x,y,z,t)u{x,y,z,t) dt 

co 1 J 


(38) 


where u(x,y, z,t) is the acoustic velocity. Assuming 
that the processes are ergodic in correlation , this tem- 
poral average may be replaced by an expected value 

I(x,y,z) = {p(x,y,z,t)u(x,y,z,t)) 

= E\p(x,y,z,t)u(x,y,z,t)] 

In general, the relationship between the acoustic pres- 
sure and velocity is not known. However, in the far- 
held, the radial velocity is dominant and the propagat- 
ing waves appear planar. Consequently, the relation- 
ship between the acoustic pressure and velocity can 
be approximated by 20,21 


u(R, 9, <j>, t) « u r (R, 9, <j>, t) = — — 2 — (40) 

pc 

where the spatial description has been changed to 
spherical coordinates for consistency with equation 
37. Thus, the far-held transmitted acoustic intensity 
becomes 


= —E[Pt(R’ 8,4>,t)] = ^ (41) 

pc pc 

where the mean square pressure has the form given in 
equation 37. 

The acoustic power incident upon a surface is 
defined as the integral of the normal acoustic intensity 
over the surface.; 


n s = j I(R,9,<f>) ■ n dS (42) 

s 

where n is the surface unit normal. Therefore, the 
total acoustic power transmitted by the plate can be 
calculated by integrating the acoustic intensity over a 
far-held hemisphere as follows. 


2ir f/2 

n t = / / It(R, 9, <j>) R 2 sin 9 d9 d<j> (43) 

<t>=o e=o 


where the transmitted intensity has the form given in 
equation 41. This integral must be evaluated numeri- 
cally, which can be accomplished via Gauss-Legendre 
numerical integration . 

Recall that the incident acoustic pressure consists 
of obliquely incident plane waves. The incident acous- 
tic intensity is then, from the above discussion, 


h = 


pc 


(44) 
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and the incident power is simply the normal acoustic 
intensity multiplied by the area of the plate: 


n,- = 


pjab cos 9i 
pc 


(45) 


Finally, the transmission loss (TL) from the incident 
to the receiving acoustic field is defined as 


TL = 10 log ( — ) 


(46) 


Results and Discussion 

The results presented in this paper were generated 
using the following material properties: 

Aluminum 


E = 68.9 GPa (10 Msi) 

G = 25.8 GPa (3.75 Msi) 
v = 0.33 

p = 2751.4 Ivg/rn 3 (0.2588xl(L 3 lb-s 2 /in. 4 ) 
Graphite-Epoxy 


Ej = 155 GPa (22.5 Msi) 

E 2 = 8.07 GPa (1.17 Msi) 

G 12 = 4.55 GPa (0.66 Msi) 

Vi 2 — 0.22 

p = 1550.07 Ivg/rn 3 (0.1458x1(7 3 lb-s 2 /in. 4 ) 
«i = -0.07 x 10 45 /°C (-0.04x10 4 7°F) 
a 2 = 30.1 x 10 45 /°C (16.7 x 10 43 /°F) 

SMA (Nitinol) 


T s = 37.78 °C (100 °F) 

T f = 62.78 (145) 

E From Fig. 9 of Ref. 3 

u* From Fig. 8 of Ref. 3 


G = 24.86 MPa (3.604 Msi) for T < %. 

= 25.6 MPa (3.712 Msi) for T > T f 

v = 0.3 

p = 6450 Ivg/rn 3 (0.6067xl0 4i lb-s 2 /in. 4 ) 
a = 10.26 xl0 4 7°C (5.7xl0 4 V°F) 

A simply supported aluminum panel (15^6x0.032 
in.) was analyzed for comparison to results presented 
by Roussos 1 . A comparison of the transmission loss 
(TL) predicted using each method for an incident 
plane wave at ${=60° and </,;=0° is shown in figure 3. 
It can be seen that the present method agrees very well 
with the classical approach for excitation frequencies 
up to 2 kHz. The deviation above 2 kHz is attributable 
to the modal resolution of the finite element model (el- 
ement size 1.07x1 in.). The radiation predictions are 
only accurate to within the ability of the finite ele- 
ment model to capture the panel vibration response, 
as is the case with any other approach. 

Transmission loss curves were also generated 
for simply supported (0,45,-45,90 ) s composite panels 
(15x6x0.032 in.) with and without shape memory 
alloy wire reinforcement. The transmission loss as a 
function of harmonic excitation frequency is shown in 
figure 4 for two panels at ambient temperature (70° F); 
a graphite-epoxy panel without SMA reinforcement, 
and a hybrid graphite-epoxy panel with 10% SMA vol- 
ume fraction and 5% initial strain in all layers (SMA 
inactivated). The shift in the first transmission loss 
dip at the fundamental frequency, and the moderate 
TL increase for most of the frequency range above the 
fundamental frequency, is attributable to the increased 
mass of the SMA reinforced panel. 

The transmission loss for the same hybrid com- 
posite panel at 150°F and 300°F (i.e. , SMA activated) 
is shown in comparison to the panel with inactive SMA 
in figure 5. It can be seen that at the elevated tem- 
peratures the increased stiffness due to tensile inplane 
recovery forces has greatly increased the fundamen- 
tal frequency, and thus the fundamental-frequency TL 
dip. Consequently, there is a significant increase in the 
transmission loss at low frequencies (below 100 Hz). 
However, at higher frequencies the TL is nearly un- 
changed. This is due to the fact that the TL is mass 
controlled at high frequencies. The reduction in the 
low frequency TL for the hybrid panel at 300° F, rel- 
ative to that for the 150°F panel, is caused by coun- 
teracting compressive inplane forces from thermal ex- 
pansion effects. 

It appears from the above discussion that there 
are two ways of improving the TL for SMA hybrid 
composite panels: optimizing the use of the SMA 
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fiber reinforcement for maximum stiffening effect with 
minimal mass addition, and reducing the radiation 
efficiency of the low order modes (particularly the 
first). These observations are in agreement with re- 
sults presented by Liang et. al. 14 and Anders et. al. 15 , 
who have initiated efforts addressing the radiation ef- 
ficiency problem. The increased flexibility of the finite 
element structural basis for the present analysis may 
prove beneficial in future studies aimed at exploring 
methods of maximizing the SMA recovery effect and 
altering the vibration mode shapes so that their radi- 
ation efficiency is reduced. 

Conclusions 

A method has been developed for predicting the 
vibration response and resulting far-held acoustic ra- 
diation for composite panels subject to elevated tem- 
peratures and an obliquely incident acoustic excita- 
tion consisting of plane waves of arbitrary temporal 
characteristics. The panel was modeled by a finite el- 
ement analysis and the radiated held was predicted 
using Rayleigh’s integral. The approach has the abil- 
ity to include other effects such as shape memory alloy 
fiber reinforcement, large deflection thermal post.buck- 
ling, non-symmetric SMA distribution or lamination, 
etc. 

Analytical results were generated for an aluminum 
panel subject to harmonic acoustic excitation and 
compared with a classical analysis for validation of 
the approach. The two approaches agreed very well. 
Results were also presented for graphite-epoxy panels 
with and without shape memory alloy fiber reinforce- 
ment. It. was found that the main effect of shape mem- 
ory alloy reinforcement on the radiated acoustic power 
is a shift in the first, transmission loss dip, due to a shift 
in the panel fundamental frequency. It was observed 
that further improvements in the panel transmission 
loss could be realized with more effective use of the 
SMA reinforcement to maximize the fundamental fre- 
quency shift and reduce the radiation efficiency of the 
fundamental mode. 
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Figure 1: Geometry on incident side of panel (after 
Roussos 1 ). 



Figure 2: Geometry on radiating side of panel (after 
Roussos 1 ). 



Figure 3: Transmission loss for a 15x6x0.032 in. alu- 
minum panel with t?* = 6 0 0 and (/>{= 0°. 
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Transmission Loss, dB 



Frequency, Hz Frequency, Hz 

Figure 4: Transmission loss for a 15x6x0.032 in. Figure 5: Transmission loss for a 15x6x0.032 in. 

(0/45/-45/90) s composite panel with 6{= 60° and (0/45/-45/90) s composite panel with 6{= 60° and 

4>i=0°. <j>i=0° • 
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